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The observed spectrum of Galactic cosmic rays has several exciting features such as the rise in 
the positron fraction above ~ 10 GeV of energy and the spectral hardening of protons and helium at 
> 300 GeV/nucleon of energy. The ATIC-2 experiment has recently reported an unexpected spectral 
upturn in the elemental ratios involving iron, such as the C/Fe or 0/Fe ratios, at energy >50 GeV 
per nucleon. It is recognized that the observed positron excess can be explained by pion production 
processes during diffusive shock acceleration of cosmic-ray hadrons in nearby sources. Recently, 
it was suggested that a scenario with nearby source dominating the GeV-TeV spectrum may be 
connected with the change of slope observed in protons and nuclei, which would be interpreted as 
a flux transition between the local component and the large-scale distribution of Galactic sources. 
Here I show that, under a two-component scenario with nearby source, the shape of the spectral 
transition is expected to be slightly different for heavy nuclei, such as iron, because their propagation 
range is spatially limited by inelastic collisions with the interstellar matter. This enables a prediction 
for the primary/primary ratios between light and heavy nuclei. From this effect, a spectral upturn 
is predicted in the C/Fe and 0/Fe ratios in good accordance with the ATIC-2 data. 

PACS numbers: 98.70.Sa, 96.50.sb, 25.40.Sc 


Introduction — The observed cosmic-ray (CR) spec¬ 
trum has several unexplained features such as the 10- 
200 GeV rise of the positron fraction e+/(e“ -Fe+) [1, 2] 
and the spectral hardening of proton and helium above 
~ 300 GeV/nucleon of energy [3-5] , that are now be¬ 
ing investigated with high precision by the AMS exper¬ 
iment [6, 7]. Recently, a puzzling spectral upturn has 
been reported by the ATIC-2 experiment for nuclear ra¬ 
tios involving iron, such as the C/Fe or 0/Fe ratios at 
^ 50 GeV/nucleon of energy [8, 9]. In the traditional 
descriptions, primary CRs such as electrons, protons. 
He, C-N-0, or Fe nuclei are injected in the interstellar 
medium (ISM) by a continuous distribution of Galac¬ 
tic sources, after being accelerated to power-law spectra 
^ E~'', with V « 2-2.4, up to PeV energies. Their spec¬ 
trum is steepened by diffusive propagation in the Galactic 
halo (typical half-size L ^3-IOkpc) with diffusion coef¬ 
ficient K oc , where <5 ^ 0.3-0.7. Interactions of CRs 
with this gas of the Galactic disk (half-size h ~ 100 pc) 
give rise to secondary particles such as or Li-Be-B 
nuclei, that are expected to be A'^-times steeper than 
primary CRs. The several models based on this picture 
agree in predicting smooth power-law spectra for primary 
nuclei, almost energy-independent primary/primary ra¬ 
tios, and a positron fraction decreasing steadily as ~ E~^ 
[10]. The recently observed features in CR leptons, pro¬ 
tons, and heavier nuclei are clearly at tension with these 
predictions. The observed positron excess requires an ad¬ 
ditional leptonic component that may come from nearby 
exotic sources, such as dark-matter particles annihilation, 
or known sources, such as pulsars or old supernova rem¬ 
nants (SNRs) [2]. In the old SNR scenario, the excess 
is produced by interactions of CR protons undergoing 
acceleration in proximity of the shock waves [11]. The 
production and subsequent re-acceleration gives rise 
to a SNR component which is harder than that of pri¬ 


mary protons or electrons, E~’^, and may explain the 
AMS data [12]. Other secondary species, such as Li-Be-B 
nuclei or antiprotons, are also expected to be produced 
in a similar way. Assuming that the observed CR flux 
is entirely provided by this type of sources, this mecha¬ 
nism predicts a rise of the B/C ratio at ^ 100 GeV per 
nucleon [13]. However the measured B/C ratio does not 
show such a feature [14]. 

In Tomassetti & Donato [15], we have shown that the 
old SNR scenario is incomplete in order to account for 
the observations of CR hadronic spectra at TeV-PeV en¬ 
ergies, because these energies can be only attained with a 
magnetic field amplification mechanism which, in turn, is 
not compatible with secondary production at the shock 
[2, 16]. Besides, the spectral hardening of CR proton 
and helium suggests that different types of sources may 
contribute to their flux [17]. In our two-component sce¬ 
nario, the total CR flux is described by a nearby source 
component in the ~ GeV-TeV region, arising from an 
old SNR, and by a Galactic ensemble SNR component 
(jfi , arising from younger sources with amplified mag¬ 
netic fields, in the ~TeV-PeV region. A key consider¬ 
ation is that, due to Compton and synchrotron losses, 
the propagation length is limited within a typical 
distance ~ V t'^^K oc with cooling time 

~ 300 X A“^MyrGeV“^. A nearby source (within 
a few 100 pc) seems, therefore, necessary to explain the 
GeV-TeV flux [18]. In contrast, CR protons and light 
nuclei do not experience radiative losses, so that their lo¬ 
cal flux may arise from the contribution of a larger popu¬ 
lation of Galactic sources [19]. As shown, such a scenario 
may account both for the rise in the positron fraction 
and for the decreasing of the B/C ratio. The new AMS 
proton data, shown in Fig. 1, are also well consistent with 
a smooth flux transition as described by model. 
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FIG. 1. Energy spectrum of CR protons multiplied by . The 
solid lines indicate the model calculations. The contribution arising 
from the nearby SNR (short-dashed lines) and from the Galactic 
SNR ensemble (long-dashed lines) are shown. The data are from 
AMS [6], PAMELA [3], ATIC-2 [20], and CREAM [4]. 


In this paper I show that, under such a scenario with 
nearby source, the shape of the spectral transition be¬ 
tween the two compontents has a characteristic signature 
in the spectrum of heavy nuclei, due to a combination of 
propagation and spallation effects. In fact, the propaga¬ 
tion range of heavy nuclei like iron is spatially limited 
by inelastic collisions with the ISM nuclei, which may 
prevent CRs injected from distant sources to reach the 
Solar System. To study this effect, I make use of an 
effective calculation scheme, based on the propagation 
scale length, that enables a prediction for the ratios be¬ 
tween light and heavy primary nuclei such as C/Fe and 
0/Fe. For these ratios, a remarkable spectral upturn is 
predicted at ~ 50 GeV of energy. 

Calculations — In conventional calculations based on 
the diffusion approximation, as long as all sources have 
the same spectral properties, the model predictions for 
the spectra of CR nuclei at Earth are known to be 
only barely sensitive to the exact distribution of Galac¬ 
tic sources [21]. But in the case of distinct classes of 
sources characterized by different properties, it becomes 
important to account for the SNR spatial distribution 
[19]. This is indeed the case for the propagation of heavy 
nuclei in our two-component scenario, where the total ob¬ 
served flux arise from the superposition of two classes of 
SNRs, and S'®, that inject CRs in the ISM with differ¬ 
ent spectral shape. In particular, the Galactic ensemble 
component, S®, reflects the contribution of a large-scale 
SNR population that extend to several kpc of distance. 
Thus, the fraction of these SNRs effectively contributing 
to the local observed flux depends on the propagation 
properties of the considered element. To first approxima¬ 
tion, the typical propagation scale distance of CR nuclei, 
A^P = VRTr'^P, can be estimated by using a spallation 
time-scale t®p = 7 ^^- Here, to effectively account that 
CRs interact only where they cross the disk, the matter 
density is considered as diluted in the propagation region 
by the h/L ratio [22]. In contrast to leptons, the function 
A®P for CR nuclei increases with energy and decreases 
with the mass. Roughly, the interaction cross sections 




FIG. 2. Top: Distribution function (solid line) and cumulative 
function (dashed line) of Galactic SNRs as function of the distance 
d from the solar system. Bottom: Fraction of Galactic SNRs con¬ 
tributing to the CR flux of O (solid line) and Fe (dashed line). 


increase with the projectile mass as (t®p oc [23], giv¬ 
ing A’*P(E) oc . This trend illustrates that, for 

the CR spectrum detected at Earth, heavier nuclei must 
come from sources located in nearer regions. Clearly, this 
reduces the fraction of Galactic sources that effectively 
contribute observed flux at Earth. In order to estimate 
this fraction for the relevant CR species, one has to model 
a realistic distribution for the Galactic SNRs as function 
of the distance to the Earth. For this purpose, I follow 
closely the effective approach of Ahlers, et al. [24], where 
the spatial distribution of SNRs is determined by a toy 
Monte Carlo generation of randomly distributed sources 
drawn from a probability density function. The input 
function assumes a four-armed Galactic structure, which 
has been confirmed by a recent analysis [25]. The cal¬ 
culation provides the SNR distribution function in terms 
of Earth centered coordinates and integrated over the 
polar coordinate. The normalized probability density of 
SNRs, P{d), is shown in Fig. 2 (top) as a function of the 
distance d from the Earth. As seen, the contribution is 
probabilistically suppressed within ~ 1.5 kpc due to the 
interarm position of the Solar System; but this does not 
prevent the factual occurrence of one (or few) anomalous 
SNR event in the Solar neighborhood that would man¬ 
ifest itself as a distinctive component of the local CR 
spectrum. In the figure, it is also shown the cumula¬ 
tive fraction of SNRs falling within a certain distance d, 
F{d) = P{l)dl. From this information, the fraction of 
Galactic SNRs contributing to the CR flux detected for 
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a j-type element is estimated as F{Xj’), where can 
be expressed as function of energy or rigidity. 

The model setup follows closely our earlier work 
[15, 26]. I briefly outline the key parameters. For the 
local SNR component, the magnetic field is B =1 
and the upstream fluid speed is ui =5xl0^cms“^. 
The SNR age is t®"*' = 50 kyr. Its maximum rigidity 
is =1TV. A damping factor kb = 16 is used 

to enhance the (otherwise) Bohm-like diffusivity at the 
shock. These properties are typical for SNRs at their 
late evolutionary stages. For the large-scale population 
of the Galactic ensemble, represented by younger SNRs 
with strong shocks and amplified magnetic fields, typi¬ 
cal parameters are ui ^ 10®cms“^, B/kb ^ 100^G, and 
~ 5 PV. No secondary production occurs in SNRs 
with these properties, and their spectrum, ~ is 

independent on the exact values of the environmental pa¬ 
rameters. The spectral indices are taken as v =2.2 and 
2.1 for Z = 1 and Z > 1, respectively, while the spectra 
from the old SNR component are softer by 0.5 for all el¬ 
ements. The spectral indices of the SNR ensemble agree 
with the basic DSA predictions and with y-ray observa¬ 
tions of young SNRs [5]. On the contrary, softer spectra 
may arise from weak shocks or from environmental effects 
such as interaction between shock and dense gas or turbu¬ 
lence damping, which may well be the case in old SNRs. 
The elemental dependence of the spectral indices is a 
known feature of the GR spectrum, possibly ascribed to 
a M/Z-dependent injection efficiency in SNR shocks [27]. 
The source abundances and the cross sections for destruc¬ 
tion/production processes are those adopted from previ¬ 
ous studies [26, 28]. The diffusion coefficient is taken as a 
universal function in rigidity, K{R) = I3Kq{R/Rq)^, with 
Kq/L =0.1/5kpcMyr“^ and index i5 = 1/2, as expected 
from an Iroshnikov-Kraichnan turbulence spectrum and 
tested to the new B/C data from PAMELA. The ISM is 
assumed to be composed by 90% H and 10% He, with 
surface density 2h x n'®™=200pc x lcm“^. The solar 
modulation is described under the force-field approxima¬ 
tion [29]. The proton spectrum is shown in Fig. 1. using 
a modulation potential = 800 MV to describe the new 
AMS data. The spectrum is described a superposition of 
two source components cfip = (fp 4>p , shown as dashed 
lines, that amounts to 85 % for the nearby SNR and 15 % 
for the ensemble, at 1 GeV/n. 

To account for the propagation/spallation effect de¬ 
scribed above, the source term for the Galactic SNR 
component, S^{E), is replaced by the effective source 
term Sf = Fj{E) x Sf{E) for all j—types nuclei, where 
Ej{E) = E{XY{E)), is the fraction of Galactic SNRs 
contributing the nuclear species j-th. In the high-energy 
limit one has —>■ for all species, but lighter CR 

nuclei experience a more rapid convergence than heav¬ 
ier nuclei, due to interactions. Typical cross sections 
for collisions with the ISM are ~40mb for protons, 
^300 mb for C-N-0, and ~ 900 mb for Fe. The function 
E{R) is shown Fig. 2 (bottom) as function of rigidity for 
C, 0 and Fe. 
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FIG. 3. Energy spectra of O and Fe multiplied by , and 
nuclear ratios C/Fe and 0/Fe as function of kinetic energy per 
nucleon. The solid lines indicate the model calculations. The con¬ 
tributions arising from the nearby source (short-dashed lines) and 
from the Galactic ensemble (long-dashed lines) are shown. The 
data are from ATIC-2 [8,9], CREAM [30] , CRN [31] , and TRACER 
[32, 33]. The TRACER data on nuclear ratios are those obtained 
in Ref. [9]. Standard model predictions are also shown for the C/Fe 
and 0/Fe ratios (dotted lines). 


Results — The model predictions are shown in Fig. 3 
for the 0/Fe and C/Fe ratios and for the spectra of Fe and 
0. The two flux components are shown as dashed lines. 
The C spectrum, not shown, is very similar to the 0 spec¬ 
trum. Owing spallation, the Fe spectrum of the Galactic 
SNR population is slightly reshaped due to a “missing 
flux” from distant SNRs that do not contributed to its 
total flux at Earth. This effect is maximized in the pri¬ 
mary/primary ratios C/Fe and 0/Fe that, as seen in the 
figure, experience a remarkable spectral upturn above a 
few tenths of GeV/nucleon energies. The ATIC-2 data 
are described very well by the model. These features are 
also present in the TRACER data, as shown in Panov et 
al. [9], but not in the CREAM data. For the two ratios 
the upturn is similar and, in fact, the C/0 ratio is fea¬ 
tureless. At E ~ 1-10 GeV/nucleon, the flux is entirely 
dominated by the nearby component. At energies above 
~ TeV/nucleon, the spallation effect vanishes and the ra¬ 
tios become asymptotically representative of the spec¬ 
tral properties of the Galactic ensemble. For reference, 
the C/Fe and 0/Fe ratios arising from standard calcula- 
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tions, i.e., using one class of sources, are plotted as dot¬ 
ted lines. As discussed, conventional models are unable 
to describe any spectral change on these ratios. From 
this mechanism, similar features are expected for other 
ratios such as Ne/Fe, Mg/Fe, or Ar/Fe. Interestingly, an 
upturn in the Ar/Fe and Ca/Fe ratios was observed by 
HE AO 3 [34]. However, these elements require more re¬ 
fined elaborations due to the presence of secondary com¬ 
ponents in their flux. Also the effective approach used 
here, only suitable for primary/primary nuclear ratios, 
suffers from several limitations. For instance, the dilu¬ 
tion factor hjL used to estimated the average interaction 
rate is probably a too crude simplification that leaves an 
uncertainties on the absolute scale of A®p. In fact, since 
the observed CRs are both injected/detected from/in the 
disk, it is improbable that during their past history they 
have spent much time in the halo. Furthermore, due 
to possible inhomogeneous diffusion or convection pro¬ 
cesses, the CR transport may be either more confined in 
the disk or swept out in the halo, respectively [19, 22, 28]. 
A rigorous treatment of the problem has to account for 
all these unknowns which, however, require the use of 
better quality CR data. In particular, the Fe spectrum 
deserves more experimental investigation. Fortunately, 
we are in a proficient era for CR physics. Ongoing and 
planned space experiments such as AMS, ISS-CREAM 
[35], or CALET [36] will measure these elements over a 
large energy range. 

On the nearby source — In this interpretation of the 
C/Fe and 0/Fe ratios, the presence of a source placed 
near the Solar System is a key ingredient. From the 
model presented here, such a source is identified as a local 
SNR {d ~ few 100 pc) with low magnetization {B ^ /iC), 
slow shock speed (ui ^5x lO’^cms"^), and a high gas 
density in comparison to that of the ISM n ~lcm“^. 
The accelerated spectra are rather steep {i/ ~ 2.7) and 
limited to a maximum rigidity ^ TV. These prop¬ 
erties are appropriate for old SNRs of type la. Rem¬ 
nants of this type may be not be detectable in y-rays 
any longer, unless the rate of p-p collisions is enhanced 
by the presence of denser media such as molecular clouds. 
In this case the emission might be sufficiently high to be 
detected by the Fermi-LAT observatory. Possible exam¬ 
ples are SNRs W44, W28, W51C or IC-443 [37-40]. Indi¬ 
cations of nearby sources in the CR spectrum are found 
in several recent studies [17, 27, 41-43], and noticed from 
independent studies in connection with the local bubble 
[41, 44, 45]. Concerning the secondary production mech¬ 
anism, antiprotons are also emitted from such a SNR but. 


similarly to the case of Li-Be-B nuclei [15], no striking 
signatures are expected after accounting for both source 
components. In fact, the p/p ratio “excess” predicted 
in related studies [12, 46, 47] arises from one-component 
scenarios. The diffuse y-ray emission can be used to test 
models involving nearby sources. Roughly, one may ex¬ 
pect a diffuse spectrum which is harder than that pre¬ 
dicted by standard models, at least in the Galactic plane 
where the emission is dominated by the production 
from p-p collisions. However, the appearance of individ¬ 
ual sources in the CR spectrum demands a different cal¬ 
culation scheme, possibly beyond the usual steady-state 
description [48-50]. This description might be, in fact, 
an oversimplification of reality which does not reflect the 
stochastic nature of SNR events and their influence on 
the surrounding CR flux. Interestingly, indications of 
CR flux variations in the Galaxy are found by recent 
Fermi-hPPT observations of the diffuse y-ray emission 
[51], and were also noted in previous studies [52, 53]. 
The Fermz-LAT data show harder y-ray spectra in the 
inner Galaxy which are at odds with standard calcula¬ 
tions. While these observations might be explained in 
terms of SNR properties [51, 53], a complete calculation 
-possibly accounting for their space/time discreteness- 
has never been attempted for y-rays. 

Conclusions — This work is motivated by the search of 
a comprehensive, agreeable model for Galactic CRs that 
is able to account for the several puzzling features re¬ 
cently observed in their spectrum. Without the presence 
of nearby sources, it is difficult to interpret the ATIC-2 
data in the context of known models of CR propagation. 
The only alternative interpretation of the spectral up¬ 
turn is the one proposed by the ATIC-2 Collaboration, 
the closed Galaxy model with bubbles [9]. In my opinion, 
their model may represent a viable solution of the puz¬ 
zle, but it suffers from some problems. In particular, it 
predicts a too weak rise of the elemental ratios, while re¬ 
quiring too steep source spectra and a too flat B/C ratio. 
On the contrary, the interpretation presented here ac¬ 
counts well for the basic observations on primary spectra 
and secondary/primary ratios. Despite the approximate 
calculation method employed in this work, it is remark¬ 
able that the spectral upturn arises under a scenario that 
is able to simultaneously account for important features 
in the CR spectrum, namely, the rise in the positron frac¬ 
tion and the spectral hardening of proton and nuclei. 
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